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Abstract: We extend our studies of holographic entropy inequalities to gapped phases of 
matter. For any number of regions, we determine the linear entropy inequalities satisfied by 
systems in which the entanglement entropy satisfies an exact area law. In particular, we find 
that all holographic entropy inequalities are valid in such systems. In gapped systems with 
topological order, the “cyclic inequalities” derived recently for the holographic entanglement 
entropy generalize the Kitaev-Preskill formula for the topological entanglement entropy. 
Finally, we propose a candidate linear inequality for general 4-party quantum states. 
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1 Introduction 

In recent years, the study of quantum entanglement and quantum information in general 
has produced a myriad of applications in high energy physics and condensed matter physics. 
A key tool for the quantihcation of entanglement is, in particular, the entanglement entropy. 
For general quantum systems, the von Neumann entanglement entropies of subsystems are 
known to obey subadditivity, the Araki-Lieb inequalities, weak monotonicity, and strong 
subadditivity, respectively: 

Subadditivity: 

S{A) + S{B)-S{AB)>0. (1.1) 

Araki-Lieb: 

S{C) + S{ABC) - S{AB) > 0. (1.2) 

Weak Monotonicity: 

S{AB) + S{AC)-S{B)-S{C)>0. (1.3) 

Strong Subadditivity: 

S{AB) + S{AC)-S{A)-S{ABC)>0. (1.4) 

Such inequalities are important, as they constrain the phase space of entanglement 
in quantum systems and can in turn be translated into other physical quantities. In par¬ 
ticular, in condensed matter physics there exists a conjectured relationship [1, 2] between 
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the existence of a gap in a system and whether or not the entanglement entropy in that 
system obeys an area law [3]. It therefore seems a fruitful direction to explore and better 
characterize the properties of the entanglement entropy in quantum-mechanical systems. 

It should be noted, however, that entropy inequalites for general quantum states are 
relatively rare; indeed, (1.1)-(1.4) are the only unconditional entropy inequalities known to 
date. Luckily, there exist classes of quantum systems for which the entanglement entropy 
is easier to characterize. 

1.1 Entropy inequalities from holography 

In holography, it has been shown that entanglement entropies of regions on the boundary 
are equal to the areas of the minimal surfaces subtending the boundary region, or in terms 
of the celebrated Ryu-Takayanagi formula [4, 5]: 

Area 

The Ryu-Takayanagi formula gives us a powerful new tool for computing entanglement 
entropies in regimes where such calculations are usually intractible. In higher dimensions, 
for example, it turns what would be a difficult (if not impossible) conformal field theory 
calculation into a straightforward minimization of area in a classical metric. 

Interestingly, these holographic entenglement entropies obey a larger set of inequalities 
than those obeyed by the generic quantum mechanical systems. In [6], it was discovered 
that, indeed, there is a new entanglement entropy inequality which is true for all systems 
with semi-classical holographic duals, i.e., the conditional mutual information is monoga¬ 
mous, or 

S{AB) + S{AC) + S{BC) > S{A) + S{B) + 3(0) + S{ABC). (1.6) 

This was done using a method known as inclusion/exclusion, in which minimal surfaces 
corresponding to the (positively signed) entropic terms on the left hand side are reparti¬ 
tioned into non-minimal surfaces corresponding to the terms on the right hand side. As 
non-minimal surfaces have more area than minimal surfaces, if such a partitioning can be 
done, then the inequality is true. A more detailed description of the methodology is avail¬ 
able in [6]. Recently, this has also been generalized to higher numbers of regions in [7] by 
converting the geometric procedure described above to a combinatoric set of contraction 
mappings from points on a hypercube, which corresponds to the left hand side entangle¬ 
ment entropies, to points on another hypercube, which corresponds to the right hand side 
entropies. This new method has yielded a new, infinite family of inequalities that has 
been proven for holographic systems. These so-called “cyclic” inequalities for n = 2k + 1 
subsystems take the form 

k—1 k 

CYC = E I{Ai...Ai : : Afc_|_j+2...A2fc+i) — ^/(Ai...Aj : : A^+fc+i) 

1=1 j=i 

n 

= S{Ai\Ai^i... Aj_|_fc) — S'(Ai... An) > 0, 
i=l 
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where I{A : B : C) = S{A) + S{B) + S{C) - S{AB) - S{AC) - S{BC) + S{ABC) and 
5(X|y) = S{XY) — S{Y) is the conditional entropy. This new generalization has also led 
to the discovery of several further holographic entropy inequalities [7]. 

Holographic systems are not the only class of quantum systems that obey a more 
restrictive set of entanglement entropy inequalities. The set of stabilizer states in quantum 
error correction does so, as well [8-10]. It is interesting, however, that the known stablizer 
inequalities are implied by (weaker than) the holographic inequalities, thus suggesting an 
nontrivial relationship between these two types of states. 

It is important to note that the utilization of holography in the inclusion-exclusion proof 
technique is actually quite minimal. Instead, holographic entropy inequalities are reduced to 
linear inequalities between the areas of boundaries of certain bulk regions. However, these 
inequalities are then proved for arbitrary bulk regions, not only for those that minimize 
the Ryu-Takayanagi entropy. As suggested in [7], it is therefore natural to expect that 
they hold likewise in condensed matter systems that satisfy an area law. The exploration 
of this idea, and the applicability of the resulting inequalities for other condensed matter 
systems, will be the focus of the remainder of this work. As any system with an exact area 
law entropy scaling necessarily satisfies these entropy inequalities, they may also provide 
further indication, particularly in the direction of falsihcation, as to whether a gapped 
system indeed implies area law scaling for entanglement entropy. 

Another potentially interesting relationship here is to the field of AdS/CMT [11, 12]; 
as we will see, the entanglement entropies for gapped phases of matter with an exact area 
scaling obey the constraints of general holographic systems, which is suggestive of possibly 
nontrivial holographic duals of condensed matter systems. 

1.2 Organization 

In this work we extend and characterize the realm of applicability of the holographic entan¬ 
glement entropy inequalities to condensed matter systems. The organization of the paper 
will be as follows: In section 2, we formally prove the validity of the cyclic inequalities for 
systems that obey an exact area law. In section 3, we show that these inequalities have a 
valid interpretation as the topological entanglement entropy in two spatial dimensions. In 
section 4, we fully characterize the entanglement entropy in systems satisfying an exact area 
law, and we give a minimal and complete set of entropy inequalities and equalities for any 
hxed number of regions. We comment about the analogous problem for general quantum 
system and propose a candidate inequality for four-partite quantum systems. Finally, we 
conclude in section 5. 

2 Gapped systems with trivial topological order 

Here we consider gapped systems with trivial topological order in d+1 dimensions. The 
entanglement entropy S{A) of a subsystem A, which measures the entanglement between 
A and its complement A^, is defined to be the von Neumann entropy 

S{A) = -tr PA log PA, 
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where pA is the reduced density operator for A obtained by tracing out all degrees of freedom 
outside of A in the many-body ground state of the gapped system. 

Note that for a bipartite system in a pure state, the reduced density operators obtained 
by tracing out either part have the same set of eigenvalues, hence the same von Neumann 
entropy, which can be seen via a Schmidt decomposition [13] of the pure state that we 
started with. This implies that for any division of the system into subsystems A and its 
complement A^, 

S{A) = (2.2) 


is satisfied. 

For gapped systems with trivial topological order, we assume that the entanglement 
entropy S{A) of a subsystem A scales as the area of its boundary and neglect any sub-area 
scaling for the moment. Namely, 

S(A) ~ dA. (2.3) 

Note that for any two regions A and B in d spatial dimensions that are non-overlapping 
except possibly at their boundaries, the entanglement entropy of the combined region AB 
satisfies 

S{AB) r-. d{AB) = dA + dB -2Ar\B. (2.4) 

Here A D B denotes the area of the codimension 1 hypersurface where regions A and B 
intersect. The above follows from (2.3) because for any d-dimensional regions A, B, and C, 
the triple intersection A n H n C is of measure zero. 

We claim that the cyclic inequality (1.7) for n = 2k + 1 regions is satisfied as a strict 
equality in this system, namely 

n 

Y, ... Ai+k) = S{A ,... An), (2.5) 

i=l 

where the sum is cyclic and all indices are taken modulo n. To prove (2.5), we compute: 

n 

LHS ~ ^ d{Ai ... Ai+k) - d{Ai+i ... Ai+k) 

i=l 

= I dAi — Y^ 2 n Aj 2 — Y^ ^ ^j 2 1 

iiiAsQo ji,j2&Qi 

n 

=dAi -|- • • • -|- dAn — ^ ^ ^ ^ 2 Aj n Aj, 

i=i jeQo 

RHS ~9(Ai ...An) 

—dA\ -j- • • • “t" dAn — ^ ^ 2 Aj.^ n Aj 2 , 

ii J2 6A/' 
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where Qi = {i + I,... ,i + k}, M = {1,...,n} and {ji,jh} are unordered pairs that take on 
the indicated values. Hence to prove that LHS = RHS, we need to show that 

n 

Y, ( 2 . 6 ) 

{jl,j2} *=lj6Qo 

ii J2 6A/' 

The LHS sums over n(n + l)/2 = (fc + l)(2A; + l) distinct unordered pairs of indices {ji, jh}) 
so it contains (k + 1)(2A: + 1) distinct terms. The summation on the RHS also contains 
{k + 1)(2A; + 1) terms. So in order to prove that LHS = RHS, it suffices to prove that any 
term appearing in the summation on the LHS also appears in the summation on the RHS. 
This is equivalent to proving that for any unordered pair {ji, jh}) ji,j 2 £ {1; • • • there 
exists i G {1,..., n},j G {i,..., i + k}, such that {ji, jh} = {l j}- 
We have the following two cases: 

1- J 2 G {ji, • • •, ji + k}. In this case, we just take i = ji,j = j 2 . 

2- 32 i {ji, ■■■,31+ fc}- In this case, if ji ^ {j 2 , ■■■,32 + k}, then {ji,..., ji + A:} n 

+ A:} = 0, otherwise ji + hi = j 2 + ^ 2 , for some /ii,h 2 G {!,..., A:}. 
Therefore, either ji = j 2 +h 2 —hi, or jh = ji+hi—h 2 . Since either hi—h 2 G {1,..., A:}, 
or h 2 — hi G {1,..., A:}, we are forced to conclude that either j 2 G {ji, ■ ■ ■ , 3 i + A:}, 
oi' ji G {j 2 ,..., j 2 + A;}. Both lead to contradictions. So if 32 ^ {ji ,..., ji + A:}, and 
if h i {h, ■■■,32 + A;}, then {ji, ...,ji + k}r\ {j 2 , • • •, j 2 + A:} = 0. Since each set 
contains A: + 1 distinct numbers, and if their intersection is empty, their union would 
contain 2A: + 2 distinct numbers, contradicting n = 2A: + 1. Thus we hnally arrive at 
the conclusion that j 2 ^ {ji,..., ji + A:} =i> ji G { 32 , ■ ■ ■ ,32 + A:}. In this case, we 
just take i = 32,3 = Ji- 


Combining cases 1 and 2, we conclude that for any unordered pair { 31 , 32 }, 3i, 32 G {1,..., n}, 
there exists i G {1,... ,n},j G {i,... ,i +A:}, such that {ji,j 2 } = {lj}- Hence (2.6) indeed 
holds and (2.5) is exactly satished by such systems. 

In section 4 below, we will generalize this result and identify all entropy inequalities 
and equalities that are obeyed in systems with an exact area-law scaling. We will hnd that 
any holographic entropy inequality is valid for systems with an exact area law. 

3 Topological entanglement entropy 
3.1 Construction and validity 

Here we consider general gapped systems in 2 -|- 1 dimensions. It is shown in [14-16] that 
the entanglement entropy of a region A with smooth boundary has the form 

Sa = OiL — 6 o 7 ( 3 - 1 ) 

in the limit L/o —)• 00 where a is the correlation length. Here bo denotes the number of 
connected components of dA, the boundary of region A. The topological entanglement 


- 5 - 


entropy —7 is a universal constant characterizing the topological state and a is a non- 
universal and ultraviolet divergent coefficient dependent on the short wavelength modes 
near the boundary of region A. In particular, 7 = logP captures the far-IR behaviour of 
entanglement and the total quantum dimension, P, which can be obtained from topological 
quantum held theory computations, is related to the number of superselection sectors of 
the system[14]. 

We divide the plane into 2k + 2 regions, labeled by Aq, Ai, A 2 , ■ ■ ■ ,A 2 k+i, where Aq 
labels the complement of A 1 A 2 ■ ■ ■ A 2 k+i-, i-e., ^0 = (^ 1^2 ■ ■ • ^ 2 fc+i)^- In order for the 
topological entropy 5topo dehned in (3.3) to be a topological invariant, we require the 
division of the plane to satisfy 

= 0, for all I C {0,1, 2, • • • , 2/c + 1}, such that 0 G I, and |/| > 3. (3-2) 

i&I 

In other words, there is no point on the plane that is shared by Aq and more than two other 
regions. We dehne the topological entropy Stopo for 2k + 2 regions as 

2 fc+l 

‘S'topo = ^ S{Ai\Ai^i... Ai+k) — <5(^1 • • • ^ 2 fc+i), (3.3) 

i=l 

where all indices are taken modulo {2k + 1). Note that our dehnition of the topological 
entropy reduces to the Kitaev-Preskill one [14] when k = 1 (i.e., a division of the plane into 
4 regions). Also note that our calculation in section 2 implies that for gapped systems with 
trivial topological order, 5topo = 0 , that is, the dependence of 5topo on the length of the 
boundaries cancels out. 

To see that S'topo is a topological invariant, consider deforming the boundary between 
regions labeled by the index set J C {0,1, 2,..., 2A: + 1}, i.e., points in the set 

s ^ n A- 

f6J 

We note the following properties of the entanglement entropy before proceeding to the 
main arguments: 

• Under general deformations of the plane, the change in the entanglement entropy 
of any region A, AS(A), is equal to the change in the entanglement entropy of the 
complement of A, AS(A‘^). This follows as a consequence of (2.2). 

• For any k ^ J, points of deformation (points in S) are far from A^. Therefore, we 
expect AS(Afc) = 0, provided that all regions are large compared to the correlation 
length. In the same spirit, AS(A) = 0 for any region A that is a union of such A^’s. 
It then follows if A is appended to any region B, the change in entanglement of that 
region is unaffected, namely, AS{B U A) = AS{B). 

Now we argue for the topological invariance of Stopo- Possible deformations of the 
regions are classihed into the following two cases: 
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1. 0 ^ J. In this case, 

2 fc+l 

^'S'topo — E \^S ^*S*(-A2 -(-i ... ^2 h-/c)] ... ^2/c+i) 

i=l 

2 fc+l 

= [A5((Ai^i+i... A^^r) - A5(^i+i... A+A - AS((Ai ... A2k+ir) 

i=l 

2 fc+l 

= ^ [A5(^o^i+i+fc^i+2+fc ■ • • ^i+2fc) — A5(Ai+i... Ai^k)] — AS'(^o) 
i=l 

2fc+l 2/c+l 

= ^ AS'(^i+i+fc^i+2+fc • • • Ai+2k) — ^ A5(Aj+i... Ai^k) 
i=l i=l 

2fc+l 2kAl 

= ^ AS'(^i+iAi+2 • ■ • Ai+k) — ^ A5'(Aj+i... Ai+k) = 0, 

i=l i=l 

where in the last step, we cyclically left permute the summands in the first summation 
by k steps, which leaves the summation invariant. 

2. 0 G J. In this case, by condition(3.2), we can either have | J| = 2 or | J| =3. 

(a) I J| = 2. Denote the only nonzero element in J as j. In this case, 

A5topo =^S{AjAjj^i... Aj^k) + Ai.S{Aj-iAj ... Ajj^k-i) + • • • 

+ ^SiyAj—kAj—k+\ ■ ■ ■ Aj^ ^S[AjAj^i... Aj^k—i) 

— AS{Aj_iAj ... Aj^k- 2 ) — ■ ■ ■ — AS{Aj_k+iAj_k+2 ■ ■ ■ Aj) 

— AS{AiA 2 . ■ ■ A2k+i) 

={k + l)ASiAj) - kASiA) - AS{Aj) = 0. 

(b) |J| = 3. Denote the nonzero elements in J as ji and j 2 . Moreover, since 
ii) j 2 G {1, 2, • • • ,2k + 1}, |ji — j 2 | < k, i.e., they are separated by a distance of 
at most k (note ji, j 2 are mod [2k + 1) integers). Without loss of generality, we 
assume j 2 = ji + I, for some 0 < I < k. We further write ji as j for simplicity. 
In this case, 

AS'topo =AS{Aj ... Aj^i... Aj^k) + • • • + AS'(Aj_|_;_fc ... Aj ... AjA 

+ AS[Aj^i-k-i. ■ - Aj ... Aj^i-i) + • • • + AS[Aj-kAj-k+i ■ • • Aj) 

+ AS[Aj^iAj^i^i... Aj^i^k) + • • • + A5(Aj+i... Aj^i... Aj^k+i) 

AS[Aj ... Aj^i... Ajj^k—i) ■ ■ ■ AS’(^j_|_;_fc_|_i... Aj ... Aj^i) 

^S[Ajj^i_k ■ ■ ■ Aj ... Aj^i_i) • • • AS[Aj_k-\-iAj_k-\-2 ■ • • Aj) 

— AS[Aj+iAj+i+i... Aj+i+k-i) — ■ ■ ■ — A5(Aj+i... Aj+i... AjA 
— A5(^i742 ... ^2fc+i) 

={k-l + l)AS[AjAj+i) + lAS{Aj) + lAS{Aj+i) 

-[k- l)AS{AjAj+i) - lAS{Aj) - lAS{Aj+i) - AS{AjAj+i) = 0. 
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To see that 5topo is a universal quantity, we consider the same argument in [14], where 
a smooth deformation of the local Hamiltonian during which no quantum critical points are 
encountered. Since the Hamiltonian is local, any smooth deformations of the Hamiltonian 
can be written as a sum of smooth deformations of local terms. Moreover, by utilizing 
the topological invariance of Stopo) we may deform the regions in the following ways while 
keeping S'topo invariant: 

• Stretch the boundaries of the regions so that L —)• oo, and the entanglement entropy 
of a region takes the form of (3.1). 

• Deform the boundaries of the regions so that all deformation of the Hamiltonian 
happens locally in the bulk of the regions. 

We further assume that the correlation length remains small compared to the size of the 
regions throughout the deformation. Hence any local deformations of the Hamiltonian in 
the bulk only has miniscule effects for the ground state near the boundary. As a result, the 
entanglement entropies of the deformed regions (and hence S'topo) should not be affected 
by such local deformations of the Hamiltonian. 

Thus we conclude that the topological entropy we defined in (3.3) is both a topo¬ 
logical invariant (invariant under deformations of the boundary of the regions that keep 
the topology of the regions unchanged) and a universal quantity (invariant under smooth 
deformations of the Hamiltonian during which no quantum critical points are encountered). 

For a general division of the plane into 2k+ 2 regions that satisfies condition (3.2), we 
can compute the topological entropy: 

2k+l 

Stopo — ^ ^ ■ ■ ■ ^i+fc)] &o[<9(^i-|-l ■ • • ^i+fc)]^ T bo[0(^Ai . . . A2k+1 )l}. (3.4) 

i=l 

where all indices are taken modulo {2k -|- 1), and bo[dA] denotes the zeroth Betti number 
(the number of connected components) of the boundary of a region A. 

Hence Stopo is proportional to the topological entanglement entropy — 7 , with the pro¬ 
portionality constant determined by the topology of the regions. This implies that we can 
extract the topological entanglement entropy of a 2 -|- 1 dimensional topologically ordered 
system with a mass gap by suitably divide the system into 2 fc -|- 2 regions, and compute the 
topological entropy Stopo- 

3.2 Examples 

Here we consider a few examples which elucidate some of the general constructions in 
section 3.1. Figures 1 and 2 illustrate two possible divisions of the plane into 2A:-|-2 regions 
that satisfy condition (3.2). 

There are three types of deformations to the regions for both divisions. We consider 
the change in Stopo under such deformations. 

First, consider deforming the boundary between two regions in figures 1 and 2. This 
can either be the boundary between two slices of the pie, say Ai and A 2 , or the boundary 
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between a slice of the pie and the complement of the pie, say Ai and Aq. In the former 
case, region is not involved in the deformation, so case 1 of our general analysis for the 
topological invariance of S'topo implies that AStopo = 0. In the latter case, there is one 
more region (^i) besides thnt are involved in the deformation, so case 2 a of our general 
analysis implies A5topo = 0 . 

Next, consider deforming a triple point where three regions meet. Without loss of 
generality, consider the three regions Ai, A 2 and ^o- There are two more regions {Ai and 
A 2 ) besides Ao that are involved in the deformation, so case 2 b of our general analysis 
implies that A5topo = 0 . 

Note that for the division in figure 1, one could also deform the center of the pie chart, 
which is seemingly different from other points in the plane. All 2A; + 1 regions but Aq 
are involved in the deformation. However, since Aq is not involved, case 1 of our general 
analysis still applies in this case, and A5topo = 0 . 

To compute the topological entropy S'topo for these two divisions, we apply the general 
formula (3.4). For pie-chart divisions in figures 1 and 2, 


bo[d{Ai... Ai+k)] = bQ[d{Ai+i ... Ai+k)] = 1 , 


therefore the summation in (3.4) gives zero, and one simply counts bo for the boundary of 
the union of all 2 A: -|- 1 regions, which yields —7 and —27 respectively. 



Figure 1. A pie-chart division of the plane into 2fc -|- 2 regions, labeled by Aq, Ai ,..., A 2 k+i- 


3.3 Beyond area-law scaling 

In the above sections, we have considered systems where the entanglement entropy is in 
the form of (3.1). In general, (3.1) will be supplemented with various corrections. We here 
consider a few examples and examine the behavior of the topological entropy (3.3) as a 
function of k in systems that deviates from exact area-law scaling. 

Since the correction leads to imperfect cancellation of the local contributions to en¬ 
tanglement entropy, (3.3) thus yields the topological entropy up to some local correction 
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Figure 2. A pie-chart division of the plane into 2fc -|- 2 regions, with a hole in the middle, labeled 
by Aq, Ai ,..., A 2 k+i- 



k 

Figure 3. For systems that deviate from exact area law, the general behaviors of corrections to 
Stopo in the forms of l/£ and £log£ are sketched in purple and blue respectively (color online) for 
i? = 10. The numerical values are up to some unknown constant of order e or /3. 

factors. For the following examples, we restrict ourselves to the pie-chart division in 2 -|- 1 
dimensions in figure 1 with radius R in units of the lattice size a. For the sake of simplicity, 
the 2 A; -|- 1 slices are divided evenly. 

For a generic 2 -|- 1 dimensional gapped system with non-trivial topological order, the 
entanglement entropy of some region A with perimeter L large compared to the correlation 
length is given by 

Sa = OiL — 6o7 + • • • (3-5) 

where the correction from (3.1) assumes the form of /3p/LP, for all odd integers p. 

The 0{L~'P) correction to the topological entropy (3.3) is given by 
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where x = (1 + ‘Ij'K^k + I/tt. For A: > 1, (5|AiSf^pQl/SA:) < 0 for each p. The higher k 
expressions are therefore slightly less sensitive to deviations of the local piece from perfect 
area-law scaling. 

More generally, we may consider entropy scaling Sa = oL — ^ + ef{L/a) where we 
recast other small deviations in the local piece of entanglement entropy into the form of 
f{L/a). Let i = Lja for a lattice with spacing a, we here briefly sketch the behaviors 
for corrections f{i) = log(^) and f{i) = ilog{i) [17, 18]. Note that near criticality [18], 
(.\og{t} scaling becomes dominant in the local piece of entanglement entropy. For systems 
sufficiently far from a phase transition, we can treat them as an order e correction in the 
area-law systems. 

For f{t) = log(^), 

r 1 / l\2A:+l'| 

AS„p„ = *log{—(l + -) }. (3,7) 

where x = (1 -|- 2l'K')k -|- l/vr. 

And for f{i) = ilog{i), 


AS'topo = e log I ^ 


1 




ru\-\ ^ n> 


In such cases, (5|AS'topoI/5A:) > 0 for regions with boundary much larger than the cor¬ 
relation length, which renders higher k definitions of the topological entropy Stopo more 
sensitive to local corrections to entanglement entropy. In principle, for some realistic con¬ 
densed matter system with non-trivial topological order, the generalized definition (3.3) 
offers a wider range of selection where one can choose the optimal k for purposes of study¬ 
ing both the topological entanglement entropy and local deviations from area-law scaling. 


4 All entropy inequalities for systems with an exact area law 

We will now derive the full set of constraints satisfied by the entanglement entropy in 
systems with an exact area law, S{A) ~ dA. We begin with a useful construction from 
[7] that allows us to reduce from continuous geometries to a combinatorial problem. For 
simplicity, we shall assume that the system lives on a manifold without boundary. 

Let Ai,... ,An be disjoint (apart from their boundaries) regions in a system with an 
exact area law. We introduce the purifying region A^+i as the closure of {Ai U • • • U A^)^. 
The entropy of an arbitrary composite region Aj = (Jie/ ^ ^ [n-|- 1] := {1, ..., n-|- 1} 

can then be evaluated in the following way, 

S{Ai)^dAi= dAindAj, 

where denotes the complement of I in [n + 1]. 

We now consider the undirected complete graph on n -|- 1 vertices, equipped with the 
edge weights w{i,j) = dAiddAj. Let w{I,J) = denote the total weight of 
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all edges between two disjoint subsets I and J, and 6{I) = the cut function. Then 

it follows from the above that 

S{Ai) ~ dAi = 6{I). 

Conversely, for any given undirected graph with non-negative edge weights we can always 
construct a geometry and associated regions ..., ^n+i such that dAj = 6{I) (cf. [7]). 
Therefore, proving entropy inequalities for systems with an exact area law is completely 
equivalent to proving linear inequalities for the cut function. 

To determine if a linear inequality 

J2ciS{I)>0 (4.1) 

I 

holds for the cut function in an arbitrary undirected weighted graph, we expand: 


= ^c/u;(/,r) = '^ci '^w{i,j) = '^w{i,j) ^ w{i,j) 

II I i&IJ&I^^ i^j I-.i&I,j^I I:i&I xoT jGl 


In the last step, the outer sum is over edges of the undirected graph. Since the edge weights 
w{i,j) are arbitrary non-negative numbers, this immediately implies that the inequality 
(4.1) is valid if and only if 

(4.2) 


^ ci>0 

l\i£l xor jGl 


for any edge Note that (4.2) asserts simply that the inequality (4.1) holds for the 

graph with a single edge {i,j} of edge weight 1, since the cut function in this case is given 
by 


W^J)I := 6{I) 


1 z G / xor j ^ I 
0 otherwise 


(4.3) 


But note that f^ij are precisely the entropies of a Bell pair shared between subsystems i 
and j in an {n+ l)-partite pure state. In view of our reduction from systems with an exact 
area law to graphs, we thus obtain the following result: 


Lemma 4.1. An entropy inequality '^j(z[n]^l ^i^l) > 0 valid for all systems with an 
exact area law if and only if it is valid for the entropies of Bell pairs shared between any 
two subsystems Ai and Aj of the purified (n -|- 1)-partite system. 

In section 2, we had proved that the cyclic inequalities (1.7) hold with equality for 
system that satisfy an exact area law. It follows immediately from lemma 4.1 that we can 
test the validity of an arbitrary entropy equality by verifying that they hold with equality 
when evaluated for Bell pairs. In [7], it was observed that this is the case for the cyclic 
inequalities (1.7) as well as four other holographic entropy inequalities established therein. 
It follows that all these inequalities hold with equality for systems satisfying an exact area 
law. In particular, this confirms our explicit derivation for the cyclic inequalities in section 2. 
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In principle, lemma 4.1 solves completely the problem of characterizing the entangle¬ 
ment entropy in systems with an exact area law. We will now describe the set of all possible 
entanglement entropies more concretely. For this, it is useful to observe that, for any fixed 
number of regions n, the collection of valid entropy inequalities Yli — 0 

a convex cone. This cone consists of all vectors s = («S'(^/))0^/c[n] ^ formed from 

the entanglement entropies obtained by varying over arbitrary regions Ai,... ,An and all 
systems satisfying an exact area law. Following [7, 9, 19, 20], we shall call it the area-law 
entropy cone. Like any convex cone, it can be dually described in terms of its extreme rays, 
which we obtain immediately from lemma 4.1: 

Lemma 4.2. The extreme rays of the area-law entropy cone for n regions are given by the 
entropy vectors Pij of Bell pairs shared between any two subsystems in the purified (n -|- 1)- 
partite system. 

Since the entropies of Bell pairs can be realized holographically, we may think of the 
area-law entropy cone as a degeneration of the holographic entropy cone defined in [7]. It 
is arguably the smallest entropy cone that can capture bipartite entanglement. 

The following theorem then gives a complete characterization of the entanglement en¬ 
tropy in systems with an exact area law: 

Theorem 4.3. A minimal and complete set of entropy (in)equalities for systems with an 
exact area law is given by (1) the subadditivity inequality S'(^i) -|- S{A 2 ) > S{AiA 2 ) and 
its permutations, (2) the Araki-Lieb inequality S'(^i) S{Ai... An) > S{A 2 ... An) and 
its permutations, and (3) the multivariate information equalities ~ ^ 

induced by any subset V C [n] of cardinality at least three. 

Proof. We first argue that the entropy equalities in (3) are correct and linearly independent. 
Their correctness can be verified by evaluating them on the rays fdij for i < j G [n -|- 1]: 

icv r.j£icv\{i} 

By symmetry, it suffices to consider the first sum. If i ^ V then it is zero. Otherwise, 

i:ieicv\(j} J.Jcv\(ij} k=o ^ 

by the standard identity for an alternating sum of binomial coefficients, which is applicable 
since |F\{i,j}| >3 — 2 = 1. The fact that the equalities in (3) are all linearly independent 
can easily be seen by induction on |17|. 

It follows from the above that the area-law cone is contained in a linear subspace of 
dimension 2” — 1 — Ylk =3 (fc) = ^ + ( 2 ) = show that this is indeed 

the dimension of the area-law cone. For this, it suffices to observe that {(3ij)k + i/3ij)i — 
if^ij){k,l} = 2 if {i,j} = {k,l}, and otherwise zero. This not only implies that the 
many extreme rays ftij are all linearly independent, but also that the area-law entropy cone 
is cut out by the inequalities 

S{Ak) + S{Ai) > S{Aki) 
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on the subspace defined by the multivariate information equalities (3). For I < n, these are 
just the inequalities in (1), while for I = n+1 we obtain the inequalities in (2) by using the 
relation S{Aj) = S^Ajc). It is clear from the above that the entropy (in)equalities (l)-(3) 
form a minimal set. □ 

Geometrically, the area-law entropy cone is an “orthant” of dimension , as follows 
from the proof of the theorem, where we have shown that the extreme rays are linearly 
independent. We note that the set of defining (in)equalities is in general not unique as the 
entropy cone has positive codimension for n > 2. 

4.1 Generating entropy equalities from graphs 

Our method can be easily adapted to include information about the spatial connectivity 
of the regions Ai,, An that enter an (in)equality: If we can guarantee that Ai n = 0 
then we do not need to consider the corresponding Bell pair /3jj when verifying an entropy 
(in)equality by using lemma 4.1. For a concrete example, consider the conditional mutual 
information I{A : C\B) = S{AB) + S{BC) — S{ABC) — S{B), which is equal to zero for 
all Bell pairs except for I3ac- If we choose A, B, and C as in the figure 4 below then 
I{A : C\B) = Q for systems with an exact area law, since ^ H C = 0. This cancellation has 
been used in [15] to extract the topological entanglement entropy. 



Figure 4. In systems with an exact area law, the conditional information I {A : C\B) vanishes for 
this configuration of regions [15]. 

This example has the following pleasant generalization: 

Lemma 4.4. Let {V,E) be an undirected graph on the vertex set V = [n], n > 3, and 
let Ai,..., An be regions such that dAi n dAj ^ 0 only if {i,j} G E. Then we have the 
following entropy equality, 

(-1)1^1 5(J) = 0, (4.4) 

/Cy JG7ro{I) 

where 7ro(/) denotes the connected components of the induced subgraph with vertex set I. 

For a complete graph {V,E), (4.4) is precisely one of the multivariate information 
equalities proved in theorem 4-3 for arbitrary regions. 

Proof. It suffices to argue that the difference to the multivariate information vanishes given 
our assumption on the spatial connectivity of the regions A\,..., An- For this, note that. 


14 - 








for any / C f/, 


5(J) -5(/)= ^ = 

\JGtto(I) ) J^K£ttq{I) j&J,k&K 

since j and k are in different connected components so that {j, k) ^ E and therefore 
dAj n dAk = 0 by our assumption. □ 

For the graph displayed in hgure 5 below we recover the statement derived above that 
I{A : C\B) = 0 for systems with ACiC = ^. 


A 


B 


C 


Figure 5. Graph corresponding to J(A : (7|i3) = 0 for regions with A n C = 0. 

We note that the connectivity assumption in lemma 4.4 is equivalent to requiring that 
I{Ai : Aj) = 0 for all {i,j} 0 E. We may therefore think of 4.4 as a constrained entropy 
equality in the sense of [21]. Below we list all non-trivial constrained entropy equalities 
obtained by lemma 4.4 from graphs with four vertices: 

(e.g., SiAiA^A^A^) - SiA^A^A^) - S{A^A 2 A^) - SiA^A^A^) 

+ 5 (^ i ^ 4 ) + ^(^ 2 ^ 4 ) + S{A^Ai) - 5 ( 214 ) = 0 
4 

(C4 = K 2 , 2 )-. S{AiA2AsA^) - Y, 5(^,|^i+i^,+2) = 0 

i=l 

(Diamond): 5(^i^2^3^4) - 5(yliA2^3) - S{A 2 A^Ai) - 5(A3^4^i) 

- 5 (^ i ^ 2 ^ 4 ) + S{AiA 2) + S{AiAi) + S{A2A3) 

+5(^2^4) + ^(^ 3 ^ 4 ) - S{A2) - 5(^4) = 0 
{Ki = W4): S{AiA2AsA^) - S{AiA 2A3) - S{AiA2Ai) - 5 (Ai^ 3 ^ 4 ) - 5 (A 2 ^ 3 ^ 4 ) 

+5(^i^2) + SiAiAs) + SiAiAi) + 5(A2^3) + ^(^ 2 ^ 4 ) + ^(^ 3 ^ 14 ) 

-S{Ai) - 5 ( 212 ) - 5(^3) - S{A4) = 0 

The last equality is in fact unconditionally true for system with an exact area law; it is the 
fourpartite information equality from theorem 4.3. 

4.2 The search for general quantum entropy inequalities 

The simplistic method of graph combinatorics has thus far managed to reproduce the forms 
of many familiar entropic inequalities. In light of the observation, one may suspect that it 
could also be useful in generating other generic n-party quantum inequalities. In light of 
the above speculation, we attempt a simple search for 4-party linear quantum inequalities 
beyond the von Neumann entropy inequalities (1.1)-(1.4). 

To search for the 4-party non-von Neumann inequality candidates, we employ again 
the notion of an entropy cone [19, 20]. Let K = {A, B, C, D} be the labels of the 4 systems 
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A 

B 

C 

D 

E 

AB 

AC 

AD 

AE 

BC 

BD 

BE 

CD 

CE 

DE 

Family 1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

0 

Family 2 

1 

1 

1 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Family 3 

1 

1 

1 

1 

0 

2 

2 

2 

1 

2 

2 

1 

2 

1 

1 

Family 4 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Family 5 

2 

1 

1 

1 

1 

3 

3 

3 

3 

2 

2 

2 

2 

2 

2 

Family 6 

1 

1 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

2 

2 

2 

Family 7 

3 

3 

2 

2 

2 

4 

3 

3 

3 

3 

3 

3 

3 

4 

4 

Family 8 

3 

3 

3 

3 

2 

4 

4 

4 

5 

4 

4 

5 

6 

5 

5 


Table 1. Families of extremal rays of the 4-party von Neumann cone constructed using known 
quantum inequalties. For a (mixed) state with subregions A, B, C and D, the region E is the 
corresponding purifying region. 


and E be the purifying system. Given a density operator pabcd, we obtain the entropy 
vector V = (vj) G where % ^ I K and vj = S{I)p denotes the von Neumann entropy 
of the reduced density matrix p/ obtained by tracing out all but the systems in I. Consider 
the set r| C of all entropy vectors v produced by physical 4-party quantum states, we 
define the 4-party quantum entropy cone as the closure r|, which is shown to be a a convex 
cone in the entropy vector space Similarly, the von Neumann cone r 4 can be defined 
as the set of vectors that satisfy the von Neumann entropy inequalities for 4-party systems, 
i.e., positivity, strong subadditivity, and weak monotonicity. Because these inequalities 
hold for an arbitrary quantum system, it necessarily follows that C r 4 . It was proven 
that r* = r„ for n < 3. Therefore von Neumann inequalities completely characterize the 
quantum entropy cone for 3 or fewer parties. 

Due to the convexity of the cones, we know that all points inside an entropy cone can 
be written as a linear combination of its extremal rays. The entropy cone r 4 produced by 
all von Neumann entropy inequalities is known and is characterized by the extremal rays 
listed in table 1 [22]. 

It is shown in [21], however, that families 7 and 8 are not physically constructible. 
Therefore, it is suspected that T^ should be a proper subset of r 4 for n > 4 [21, 22],^ 
so that additional 4-party entropy inequalities may be needed to complete the entropy 
cone.^ By searching through the integral linear combinations of the constrained entropy 
equalities constructed from graphs as described in section 4.1 above, we have generated 
a set of inequalities that satisfy families 1 through 6 but can violate families 7 and 8 for 
certain permutations. After a cursory search, we found two such candidate inequalities (up 


^In fact, studies in the classical Shannon entropy reveal that additional entropy inequalities, such as 
Zhang-Yeung Inequality, are needed in addition to the Shannon-type entropies. 

^There is also the possibility of a characterization in terms of non-linear inequalities. However, that is 
beyond the scope of this work. 
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A 

B 

C 

D 

E 

AB 

AC 

AD 

AE 

BC 

BD 

BE 

CD 

CE 

DE 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


Table 2. A candidate extremal ray for the 4-party quantum entropy cone proposed by [22]. 


to permutations): 

S{ABD) + S{ABC) + S{CD) - S{BC) - S{AB) - S{BD) - S{AC) 
- S{AD) + S{B) + S{A) <0 

S{ABD) + S{ABC) + S{BCD) - 2S{BD) - 2S{BC) + S{CD) - S{AD) 


(4.5) 


(4.6) 


- S{AC) - S{AB) + 2S{B) + S{A) < 0 
Both inequalities, as well as the quantum analogue of the Zhang-Yeung inequality [23, 24], 

I{A : B) + I{A : CD) + 3/(C : D\A) + I{C : D\B) - 2I{C :D)>0, 

also satisfy the candidate extremal ray in table 2 for the 4-party quantum entropy cone 
r^. We note that the graph construction, as currently formulated, cannot produce the 
Zhang-Yeung inequality. 

Inequality (4.5) is known as the Ingleton inequality. It can also be written as 

I{A ; B\C) + I{A : B\D) + I{C : D) - I{A : B) > 0. 

It is known that the Ingleton inequality does not hold for general quantum states (not even 
for classical probability distributions), but that it is a valid inequality for the subclass of 
stabilizer states [ 8 , 9]. 

Inequality (4.6) on the other hand seems to be independent of the other 4-party linear 
candidates and to the best of our knowledge, has not been tested to a greater extent. 
All tests we’ve conducted so far on this inequality return the same result as the quantum 
analogue of Zhang-Yeung inequality. It will be worthwhile to generate random 5-partite 
quantum pure states and numerically check if the inequality can be violated. Note that such 
checks do not constitute a proof. However, it can be useful in finding a counterexample. 


5 Conclusion and future directions 

We here restate our findings: 

1. We have completely characterized the entropy (in)equalities obeyed by systems in 
which the entanglement entropy satisfies an exact area law. We find that such an 
entropy inequality is valid if and only if it is valid for the entropies of Bell pairs shared 
between arbitrary subsystems. In particular, all holographic entropy inequalities, such 
as the cyclic inequalities established recently in [7], are satisfied by systems with an 
exact area law. These (in) equalities may provide constraining tests to determine 
whether certain condensed matter systems satisfy an area law. 
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2. The cyclic (in) equalities in two-dimensional systems with non-trivial topological order 
can be seen as a generalization of [14] which extracts the topological entanglement 
entropy using higher number of partitions. These higher k generalizations of 5topo 
are sensitive (or insensitive) to different types of deviations from area-law scaling. 

3. A graph representation for constrained entropy equalities for systems with an ex¬ 
act area-law scaling is found. As this construction recovers a wide class of entropy 
equalities including strong subadditivity, it may be suspected that further quantum 
inequalities may also be found in the set of graph-generated equalities. Following this 
approach, we have found a candidate linear entropy inequality for general 4-party 
quantum states. 

As we have seen, the graph representation of entropies in area-law systems used in 
section 4 offers surprisingly powerful insights. In the absence of the minimization that 
appears in holography, several holographic inequalities now hold exactly as equalities for 
systems satisfying an exact area law, and we may understand the entropy cone spanned by 
the area law systems as a particular degeneration of the holographic entropy cone [7]. The 
method may also provide useful insight for the long-standing problem of finding linear 
inequalities for the entropies of general multipartite quantum states. In this regard, we 
also note that generalizations of the graph-theoretical approach is much desirable. One 
such generalization will involve constructing different graphs for a quantum state with 
holographic dual. We suspect that the geometry of AdS or its dual kinematic space can 
be effectively captured by analyzing generalized graph representations for these states. In 
particular, machineries developed in spectral graph drawing may be used to recover the 
emergent geometry for more general states. 
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